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The full analytical solution of the Schrodinger equation for the hydrogen molecular ion H% (special 
case of the quantum tree-body problem with the Coulomb interaction) is obtained first. The solution 
shows that the total wave function is a two-term function in the sense that it is a linear combination 
of the two linear independent wave functions. The two-term character of the total wave function was 
visualized in calculations of the total electron density of Hf at different internuclear separations. 

PACS numbers: 31.10.+Z 


The hydrogen molecular ion is a very impor¬ 
tant case of the quantum tree-body problem with the 
Coulomb interaction. It is well known that the variables 
of the Schrodinger equation for the H molecular ion in 
the Born-Oppenheimer approximation are separated in 
elliptic coordinates r] = (rq — r 2 )/R, £ = (?T + r 2 )/R , 
where R is the internuclear distance and n, r 2 are dis¬ 
tances between the electron and nuclei 1 and 2 , corre¬ 
spondingly PQ. In this case, the wave function for the £ 
ground state can be presented as 


* = X{Z,R)Y{t I ,R), 

where the functions X(t;,R) and Y(rj,R) satisfy the fol¬ 
lowing equations: 
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where E and A are energy and separation constant which 
have to be determined. Equation ([2]) is most important 
among these two equations because its solution is mainly 
defined the property of the molecular wave function along 
the internuclear axis. 

An analytic solution for the Eq. JT1) was given in j2], 
while an analytic solution of Eq. pF was unknown up 
to present study. However, its solution can be obtained 
by noting that Y -A e ±v ^ v when rj —► 0. Then, the 
two linear independent solutions of this equation can be 
obtained by substituting two anzatz functions into it 

Y 1 =g 1 (r t )e +VXr >, 


Y 2 = g 2 (r,)e-^\ 
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where < 71 ( 77 ) and < 72 ( 77 ) are new functions. This gives us 
two equations for < 71 ( 77 ) and < 72 ( 77 ) 
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92 = 0 . 


Their polynomial solutions can be found by applying the 
usual method used for the problem [ 21 0 ]. 

The first most important question is algebraic struc¬ 
ture of the general solution of differential Eq. Q. It is 
written as the linear combination of two linear indepen¬ 
dent solutions Yi and Y 2 : 


Y = C'i(77)e + ' /3 " 5 i(77) + C 2 {r])e-^g 2 {r]). 


This solution is symmetric with respect to the sign 
change of 77 . Therefore, it can be rewritten in the form 

Y = D 1 (abs( V ))e+^ abs ^f 1 (r,)+ 

D 2 (-abs(r,))e-^ abs ^f 2 (r]), 

which explicitly displays physical meaning of both terms. 
The first term in ([3| can be interpreted as the molecular 
wave function of cylindrical symmetry, because in inter¬ 
val [— 1 , 1 ] it has a minimum at 77 = 0 and maxima at 
77 = ± 1 , while the second one can be interpreted as the 
quasi-atomic wave function of spherical symmetry, be¬ 
cause it has a maximum at 77 = 0 . 

The algebraic structure of the wave function Y (|3j) di¬ 
rectly shows the correlation of the total wave function 
with the wave function of the united atom He + studied 
first in [5] and in many other investigations including de¬ 
tailed consideration given in |2j. Taking this into account 
we will call the total wave function a two-ternr total wave 
function because it is a linear combination of two linear 
independent wave functions. 

It is important to emphasize here that the algebraic 
structure of the total molecular wave function and its 
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FIG. 1. Electron densities in 74^" at different internuclear distances (in a.u.): (a) - 0.008, (b) - 0.012, (c) - 0.025, (d) - 2.0. 


asymptotic or correlation property in the direction of the 
united atom does not contradict the two-body cusp con¬ 
dition for many-body wave functions [7]. It is sufficient to 
note that the two-body cusp condition is correct only in 
the full coordinate space, which has a dimension of 3 N, 
where N is the number of considered particles [7]. The 
full coordinate space is a non-physical space. In three- 
dimensional physical space the total molecular wave func¬ 
tion has the asymptotic or correlation properties but no 
two-body cusp conditions. 

Numerical calculations of three-dimensional electron 
density distributions in at different internuclear dis¬ 
tances give additional confirmation that the total wave 
function is a two-term function. 

In this connection, the total electron density of the 
was calculated by the Hartree-Fock method using 
the Gaussian 94 program [5] . The calculations were per¬ 
formed with the nuclei centered 15s6p5d4/3g2/il i Gaus¬ 
sian basis set. Obtained total electron density distribu¬ 
tions in H 2 are presented on Fig. [l] 

From the fact that the total wave function is a two- 
term function follows that at short internuclear separa¬ 
tions, where the quasi-atomic wave function is dominant, 
the total electron density distribution should be in the 


form of a perturbed spheroid. Explicitly, such a form of 
the total electron density distribution is clearly observed 
on Fig. ID 1 at internuclear separation R = 0.008 a.u. The 
intermediate region begins at larger R , where the ampli¬ 
tude of the molecular wave function becomes noticeable. 
This leads to the appearance of two maxima on the total 
electron density distribution. Explicitly such a distribu¬ 
tion with two maxima at the positions of two nuclei can 
be observed on Fig. [T]d at R = 0.012 a.u. This distance 
is in excellent agreement with the investigation of the 
Bethe perturbation theory 5] in [9], where it was shown 
that noticeable deviations of this theory from the exact 
solution for begin at this distance. In the intermedi¬ 
ate region the amplitudes of both wave functions have to 
become comparable at some R and, hence, electron den¬ 
sities from both wave functions have to be observable on 
the total electron density distribution. Namely, this case 
is presented on Fig. jTJ; at R = 0.025 a.u. The electron 
density from the quasi-atomic wave function of spherical 
symmetry is located in the middle between two nuclei 
and the electron density from the molecular wave func¬ 
tion of cylindrical symmetry is well recognizable by the 
two maxima at the positions of the two nuclei. Finally, 
at large R the molecular wave function has to become 
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dominant and the distribution of total electron density 
has to have the usual form. Such a distribution of the 
total electron density can be seen on Fig. E* at R = 2.0 
a.u. 

Thus, the analytical solution of the Schrodinger equa¬ 
tion for the hydrogen molecular ion and the total 
electron density distributions presented above show that 


the total wave function is a linear combination of two- 
linear independent wave functions (a two-ternr function) 
in which the amplitudes of wave functions change with 
the variation of the internuclear separation. Near the 
equilibrium point the molecular wave function of the to¬ 
tal wave function is dominant, while at short distances 
the quasi-atomic wave function becomes dominant. 
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List of the figure captions. nuclear distances (in a.u.): (a) - 0.008, (b) - 0.012, (c) - 

Figure 1. Electron densities in H£ at different inter- 0.025, (d) - 2.0. 



